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It is conventional wisdom that staggered fermions do not feel gauge field topology. However, the response of 
staggered fermion eigenmodes to the topology of the gauge field can depend quite sensitively on the way in which 
the staggered fermion action is improved. We study this issue using a variety of improved staggered quark actions. 
We observe that the separation between the "would be" zero modes and the non-chiral modes increases with the 
level of improvement. This enables the "zero modes" to be identified unambiguously. The distribution of the 
remaining non-chiral modes is compared with the predictions of Random Matrix Theory. Satisfactory agreement 
is obtained. 



It is a long standing problem that staggered 
fermions do not feel gauge field topology on coarse 
lattices. In particular, studies of the microscopic 
eigenvalue spectrum of the staggered Dirac oper- 
ator, upon comparison with the analytic results 
predicted by Random Matrix Theory (RMT), 
show that only the trivial topological charge sec- 
tor is probed by the operator yQ. On the other 
hand, it has been suggested that lattice artifacts, 
in particular flavour changing interactions asso- 
ciated with the staggered quark action, are the 
main cause for this failure and other related topo- 
logical problems [2]. The staggered quark action 
describes four quark flavours in the continuum 
limit and the eigenvalue spectrum has a 4-fold 
degeneracy in this limit. At finite lattice spac- 
ing, the flavour changing interactions break the 
flavour symmetry and the degeneracy is lifted. 
Consequently, staggered fermions do not have ex- 
act zero modes at finite lattice spacing because 
the continuum chiral modes (if there are any) are 
scattered on the lattice. One thus expects stag- 
gered fermions to show better topological proper- 
ties if flavour changing effects can be suppressed. 

In this project, this issue is examined using a 
variety of improved staggered quark operators, 
which are designed to suppress flavour changing 
effects by smoothing out the quark-gluon inter- 
action vertex We observe that the "would 



be" zero modes are visibly decoupled from the 
non-chiral modes as the level of improvement in- 
creases. The effect of changing the lattice spacing 
on the eigenvalue spectrum is also studied. We 
again observe that separation between the "zero 
modes" and the non-chiral modes increases as one 
approaches the continuum limit. This enables the 
"zero modes" and subsequently the gauge field 
topological indices to be identified. Although the 
charge indices obtained by using different oper- 
ators do not agree on a configuration by config- 
uration basis, the charge distributions are found 
to have no significant difference. Finally, distri- 
bution of the remaining non-chiral modes is com- 
pared with the predictions of RMT. 

The unimproved staggered Dirac operator is 

D *>v = 5 ^ V ^ \ U ^ 5x +^y ~ U i(y) S *-»,y\ > C 1 ) 

where i]^(x) = (— l) a! i+-"+ a <*-i is the standard 
fermion phase. A variety of improvement schemes 
are considered here. The basic improved version 
is the 0(a 2 ) improved Asq operator which 
includes an additional 3-link Naik term and re- 
places the gauge field in (JTJ with Fat7 effective 
links (sum of the original link and the nearby 
paths, up to 7- link staples). Further improve- 
ment iterates this fattening procedure with an 
additional SU(3) unitarization step between suc- 
cessive iterations, giving the improved Asq op- 
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Figure 1. Eigenvalue spectra of different opera- 
tors at (3 = 5.85. Results are shown for 50 con- 
figurations. 



erators (UFat7)" x Asq 5 . The HYP-improved 
operators (HYP)™ [5] are constructed in the simi- 
lar fashion but only those links within the hyper- 
cube containing the original link are included in 
the fattening process. We refer the readers to the 
respective papers for further details. 

Simulations are done with the standard Wil- 
son gauge field action at three values of coupling, 
j3 = 5.85, 6.0 and 6.2. About 1000 configurations 
are generated for each /?. The lattice sizes are 
10 4 , 12 4 and 16 4 respectively so that the physical 
volumes are ~ (1.2/m) 4 in all cases. We compute 
the eigenvalues (A 2 ) and chirality (x) of the low- 
est 10 eigenstates of — D 2 for all operators listed 
above. The corresponding eigenvalues of D are 
HX. Since D 2 connects either even-even or odd- 
odd sites on the lattice, only half of the spectrum 
is computed. For comparison, the lowest 5 eigen- 
values (in each chiral sector) are also computed 
for the overlap operator [Jj on the 10 4 lattice. 

The effect of improvement on the eigenvalue 
spectrum is shown in figure^where |x| is plotted 
against A for the different operators at a fixed cou- 
pling (3 = 5.85. It is observed that the "would be" 
zero modes separate from the non-chiral modes 
with higher level of improvement. In addition, 
as the level of improvement increases, the contin- 
uum 4-fold degeneracy emerges where the scat- 
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Figure 2. Spectra at (3 = 6.0 and 6.2. 



tered eigenmodes begin to form quartets (remem- 
ber that only half of the spectrum is computed so 
the data points form doublets). Note also that at 
this coupling the Asq operator still retains large 
lattice artifacts so the "would be" zero modes 
cannot be identified. Further improvement is re- 
quired. 

We next study the dependence of the eigenvalue 
spectrum on (3. The spectra of the UFat7 x Asq 
operator at (3 = 6.0 and 6.2 are given in fig- 
ure [3 Again, the "would be" zero modes are 
well separated from the non-chiral modes and 
the continuum 4-fold degeneracy is better real- 
ized as one approaches the continuum limit. To 
quantify the separation between the "zero modes" 
and the non-chiral modes, the ratio of eigenval- 
ues between the smallest non-chiral mode and the 
largest "zero mode" is plotted in figure [3] The 
ratios increase from one order of magnitude at 
(3 — 5.85 to three orders of magnitude at f3 = 6.2 
for the improved staggered operators. Theoreti- 
cally, the ratio is infinite for the overlap operator 
because exact zero modes exist on the lattice for 
overlap fermions. It is finite here solely because of 
computational precision. Results here also show 
that improvement using fat-links or hyper-cubic 
blocking are equally efficient. 

Because of lattice artifacts, the topological 
charge indices obtained by using different oper- 
ators do not agree on a configuration by config- 
uration basis. In the present case we find that 
the charge indices determined by using different 
operators agree at about 60-70%, compared to 
28% if the values were completely random. It is 
then important to check whether the charge dis- 
tributions are also different because physical ob- 
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Figure 3. Ratio between the eigenvalues of the 
smallest non-chiral modes and the largest "zero 
mode" . 
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Figure 4. Normalized topological charge distribu- 
tions (a) UFat7 x Asq (b) HYP (c) overlap. 



servables, e.g., the susceptibility, are related to 
the average of the topological charge. The distri- 
butions obtained by three different operators are 
given in figure 21 ft can be observed that there is 
no significant difference among the results. Con- 
sequently, one would expect topological quanti- 
ties obtained by these operators to agree. 

Finally distribution of the remaining non-chiral 
modes is compared with the predictions of RMT. 
In particular, we focus on the cumulative distri- 
bution of the smallest non-chiral modes [I] . The 
results are shown in figure^] Agreement with the 
predictions of RMT is evident and is comparable 
to that obtained with overlap fermions 8 . 

In summary, we have shown that the eigenvalue 
spectrum of the staggered Dirac operator is very 
sensitive to lattice artifacts. Better continuum 
properties are obtained by using improved stag- 
gered operators. The topological charge distribu- 
tions calculated using different operators show no 
significant difference. Finally, distribution of the 
smallest non-chiral mode is compared with the 
predictions of RMT and satisfactory agreement 
is obtained. 



Figure 5. The cumulative distributions of the 
smallest non-chiral mode in different topological 
sectors. Solid lines are predictions of RMT. The 
eigenvalues are rescaled by the infinite-volume 
chiral condensate E, z = XVE, where £ is a fit- 
ting parameter here. 



Two related works with similar conclusions 
were also reported at this conference [0]. 
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